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Abstract 

We study the stability of the highest symmetric solution (Wigner-solution) of 
Dyson-Schwinger equations in chiral limit and at zero temperature. Our results con- 
firm that if the chemical potential is not very large, the QCD vacuum is in the chiral 
symmetry breaking phase and the quantum phase-transition of the chiral symmetry 
restoration is in first order. Meanwhile it seems that there is not competition between 
chiral symmetry breaking phase and color superconductivity phase since the color su- 
perconductivity phase appears only if the chemical potential is very large. Moreover, we 
propose that chiral symmetry breaking arises from the positive feedback with respect 
to the mass perturbation. 

PACS Numbers: 12.38.Lg, 11.30.Rd, ll.10.Wx, 25.75.Nq 



* Corresponding author, e-mail address: yxliu@pku.edu.cn 



1 



It has been known that the results of perturbative renormalization-group in QCD will 
encounter divergence at low energy region. This behavior indicates an important fact that 
the vacuum (Wigner-vacuum) defined within Feynman's original path-integral framework 
which possesses the highest symmetries is unstable and incorrect [1] . In addition, it should be 
emphasized that the effective action (viz. minus effective potential) based on Wigner-vacuum 
should have the same globle and modified BRST symmetry[2] as the original action[3]. Thus, 
the gaps corresponding to the running coupling coefficients of the variational one-particle- 
irreducible (1PI) vertexes, which break the global symmetry and BRST gauge symmetry, 
should definitely be zero! Nevertheless, the statement, which is correct in high energy region 
and is the key to prove gauge theory is perturbatively renormalizable|2J, often deviates 
from observations to real QCD in low energy region. Therefore, the real QCD vacuum in 
low energy channel is different from the Wigner-vacuum. Anyway, we should note that 
energy or free energy (if chemical potential /x ^ 0) of Wigner-vacuum is definitely the 
lowest in Euclidean formalism of QCD at zero temperature (temperature T = means 
a unbounded imaginal time). Thus we can't justify a phase-transition by comparing the 
energy between the symmetrical vacuum and the symmetry-broken one. Indeed, there has 
been a powerful general argument that, while ground states of a macroscopic system are 
degenerate, the true vacuum is the asymmetric one rather than the Wigner-vacuum because 
the Wigner-vacuum is extremely unstable against a perturbation (viz. symmetric system will 
undergo a dramatic transforms against an even arbitrarily small perturbation) [I]. However, 
such an idea has not yet been proved solidly. It is then interesting to study under what 
conditions the vacuums of QCD tend to be degenerate, and how to understand the phase- 
transitions in a reasonable and systematic way. It has been known that Dyson- Schwinger (D- 
S) equations approach provides a nonperturbative framework to study the vacuum properties 
of strong interaction and hadron properties in free space and to simulate the chiral symmetry 
restoration and deconfinement in the system with finite temperature and/or finite chemical 
potential[8l [5J (6J [7]. In this paper, we intend to shed some light on these questions by 
analyzing the stability of Wigner-vacuum in the framework of Dyson-Schwinger equation 
approach. 

Within Feynman's framework, one can easily derive a series of dynamical integral- 
equations (Dyson-Schwinger equations) [5] and a series of identities which come from the 
symmetries. If we combine these two series of results, we will get a theory based on Wigner- 
vacuum which gives Wigner solution of the equations for any correlations and gaps. On 
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the other hand, if we relax the restrictions of the symmetries, the Dyson-Schwinger (D-S) 
equation will be correct not only for the symmetrical vacuum but also for the asymmetrical 
vacuum (if existing!) and, besides the Winger solution, we can get a new class of solutions 
(Nambu solution) of the D-S equation for the correlations and gaps. We must note here that 
the Nambu solution is important if and only if the Wigner-vacuum is unstable under per- 
turbations, viz. vacuums are degenerate, otherwise, Nambu solution is not the real ground 
state but a dynamically stable excited state. For example, as we will see at below, the fact 
that chiral susceptibility of Wigner mass function in the chiral limit is negative indicates 
the chiral-symmetry breaking of the QCD vacuum in a natural sense, and the restoration 
of chiral symmetry is closely related to the presence of positive chiral susceptibility beyond 
some critical chemical potential. 

It is well known that D-S equation is a series unclosed equations where the equation 
for n-point Green functions depend on (n + 1) point Green functions. One should then 
truncate it with indispensable approximations before taking it to evaluate any physics quan- 
tity practically. Furthermore, in order to obtain a full and reliable understanding of QCD 
phase-diagram with respect to the medium density (or the chemical potential), we should 
investigate the stability of the Wigner-vacuum against all possible perturbations simultane- 
ously which are allowed by the truncation to D-S equation. In this work, for simplicity, we 
take the rainbow approximation and investigate the stability of the Wigner-vacuum against 
the allowed chiral-perturbation and diquark-perturbation. 

We begin with single-flavor QCD in chiral limit. The action in Euclidean space is usually 
given as 



Representing the 1PI vertexes of ipip, ipip, i/ji/j as S 1; S 2 , S3, respectively, we can write 
the dressed quark propagator without diquark component as 




(1) 



where the 7 matrixes are chosen to satisfy [7^,7^] + 



25^ and D^ = + ig A 




G = 



1 



(2) 



p • 7 - i/i7o + Si 



Meanwhile the relations between the 1PI vertexes and the full propagators can be written 



as 



G 



^, ^2,3 — G(—Yj2,z)G\. 



(3) 



1 — X^GrX^Cr 
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With these relations we can easily obtain the corresponding D-S equation 

Eias(p) = -9l J ^l^G^Ml^D^k) , (4) 

where k ~ p — q stands for the momentum transferred. Furthermore, if we are only interested 
in the condensates in chiral and color 3 diqurak channels which preserves the parity and in 
which single gluon exchange interaction is attractive [I], the 1PI vertexes in strong interaction 
matter with chemical potential fi can be expressed as 

Ei(p) = [A(p) - l]p • 7 + [C(p) - l](po - i//) 7 o + Ai(p) , (5) 

£ 2 (p) = A 2 (p)M7 5 C f , (6) 

Z 3 (p) = A 3 (p)CM l5 , (7) 

where M is a matrix in color space corresponding to the color 3 diquark channel, C is the 
charge conjugation operator, which can be given explicitly as 

M Q/3 = e laf) (8) 

Cm = 727o* T , = ^ T 7 2 7o • (9) 

It is evident that the Wigner solution is characterized by Ai )2) 3 = 0. 

Instead of exploiting Nambu solution of D-S equation, we prefer to investigate whether 
the Wigner-vacuum is stable against perturbations of the condensates A^vl/vl/, A\m M^^Cm 
and A\{Cm)M^^! . For our purpose, here we have defined the perturbed Wigner solutions 
for the gaps as 

Af r ^(p) = F 1 (p)A k 1I (10) 

AY l9ner (p) = F 2 (p)A b 2 , (11) 

AY t9ner ( P ) = F 3 (p)A». (12) 

The Fi, F 2i 3 stands for the susceptibility with respect to the chiral channel, diquark channels, 
respectively. 

From the D-S equation with infinitesimal but explicit mass term, we find that, when 
A\ ^ 0, A23 = 0, Fi satisfies equation 

f,w = 1 " 9A / Wf^Au.w + ca,( 9 )(* - w 7 "^ DAk) ' (13) 
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From the D-S equation with infinitesimal but explicit diquark term, we obtain that, when 
A2 7^ 0, A\ 3 = 0, F 2 satisfies equation 

d A q \ F 2 (q) 



(2tt) 4 lfl 2 A w (q)q- 7 + C w (q)(q - ifJ,)j 
xM 75 (7 1 7 Ad (*) . (14) 

Here Apy, CV are Wigner solutions for the A, C which satisfy the D-S equation with A 12 ,3 = 
0. After some derivation, the equations can be explicitly written as 

W[P) 3p>J (2tt) 4 ^( g )^ + C^( g )(g Q -W 

(15) 

and 



4gj / rf 4 g 2A w -(g)(po-^)g i D(M(fc) + C w (q)(po-ifi)(q -iij)D ff tl (k) 
3(p -z^J (2tt) 4 ^(g)g 3 + C&(g)(g - ^) 2 



(16) 



where D^, = D 00 - At- 



MM 



It is apparent that to solve the equation practically, one needs the effective gluon prop- 
agator as an input. We then take the effective gluon propagator as the Tubingen model [9 J 

with 

g\D(k) = 4vr 2 d— . 

By using relation 

h M ^L = --M, (18) 
2 2 3' V ; 

we could conclude that the stability of Wigner- vacuum against chiral and diquark perturba- 
tion is characterized by linear integral equations 

Flip) = 1 + 4, 2 / - A- F ^ qWk) (19) 



and 



(27r)MMg)? + CSr(g)(gb-v) 
A F 2t3 (q)D(k)[A 2 w (q)if + C 2 w {q){ql + /i 2 )] 



2 ' 3W V (2tt)* + C% 2 - /i 2 )] 2 + [2C^]* " 1 UJ 

Up to now, what we have discussed is only one flavor of quark. In fact, we can deal with 
the case of two or three flavors of quarks without any technical difficulty in chiral limit. Since 
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in chiral limit, flavor has nothing to do with dynamics, the gapped flavor channel should 
be chosen to leave the maximal unbroken symmetries [I], such as ^ i a e la ^ j@ in case 

of two flavors and ^i a e^e A%3 jsC^&jp in case of three flavors(a,/3 are color indexes and i,j 
are flavor indexes), which are the well known two- flavor color superconductivity (2CS) and 
color flavor locking (CFL) channels, respectively. However, all these considerations, which 
do not change the equations of the susceptibilities in chiral limit, are nontrivial in case of 
real mass spectrum of {u, d, s} quarks because it has strongly indicated a transition from 
the 2CS phase to the CFL phase at some relatively high chemical potential. In the present 
work, as we work in chiral limit, one-flavor should be sufficient. 

To obtain the quantitative result, we solve the D-S equation with the Tubingen model 
of the effective gluon propagator and parameters u = 0.4 GeV, d = 45.0 GeV~ 2 , with 
which the pion properties and some other low energy chiral observables have been described 
well [9 j. The obtained results of the chemical potential dependence of the Wigner solutions 
CVI/*] are illustrated in Fig. 1. The obtained chemical potential dependence of the 
Nambu dynamical mass function as the ratio A^v/A/v of the Nambu solutions is illustrated 
in Fig. 2. And the chemical potential dependence of the susceptibilities corresponding to 
the chiral quark gap and diquark gap of Wigner solution and .F2 [/■*]) are displayed in 

Fig. 3, Fig. 4, respectively. Fig. 1 shows obviously that the Wigner-solutions A(p, p) and 
C(p,p) are equal only if p = (another point shown in Fig.l is not universal for other 
momentum modes), but remarkably separated from each other once p 7^ 0. Since the A\y 
and Cw contain the crucial information of the QCD-vacuum in our current framework, we 
have calculated them as precisely as we can, rather than taking any approximation such 
as A = CfTO]. In fact, because of the obvious oscillating behavior of A w (p) and Cw{p) 
caused by chemical potential, which will be restrained by chiral gap in Nambu solutions, the 
numerical calculation for Wigner solutions is much more difficult than for Nambu solutions. 

Because the susceptibility is the response rate of the dynamical mass with respect to the 
perturbation of current quark mass, with the above obtained result we discuss the quantum 
phase-transition in QCD. At first we take, for instance, a free fermion system in chiral 
limit with F% = 1 as a starting point. Corresponding to the presence of the repulsive 
potential in chiral channel induced by mass vertex of the system, the influence of small 
and positive mass term on the system is nothing more than a gapped dispersion relation 
E = (p 2 + m 2 ) 1 / 2 . Changing the sign of the perturbative mass term just trivially results in 
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Figure 1: Chemical potential /i dependence of Winger-solutions in extremely low energy 
channel (p 2 = 0) 
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Figure 2: chemical potential dependence of the dynamical quark mass function in extremely 
low energy channel (p 2 = 0) 

a redefinition of fermion and anti-fermion, and of course the same physics. Therefore, the 
influence of small mass term on free system is really perturbative, and we can infer safely 
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Figure 3: Chemical potential \i dependence of the susceptibility F 1 in extremely low energy 
channel (p 2 = 0) 
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Figure 4: Chemical potential dependence of the susceptibility F 2 in extremely low energy 
channel (p 2 = 0) 

that the free vacuum should preserve chiral symmetry. Now, we turn to the case of QCD. 
Because QCD is an asymptotically free theory, the chiral susceptibility in high energy mode 
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is positive to maintain the system stable in high energy region. If the F\ in low energy mode 
is also positive, the QCD should be similar to the free case by means of that the vacuum 
is stable against mass perturbation. In fact, our calculation indicates that, in low chemical 
potential region, the F\ is negative in low energy channel. In this case, a small and positive 
mass term mM>^> results in an attractive (not repulsive as the free case!) potential in the 
chiral channel whose strength is characterized by |mFi|. This attractive interaction, even 
if quite weak, should induce a small (\E , \[>) condensate correspondingly. Since the induced 
small condensate can be treated as a positive mass-like vertex, which will be dressed by 
QCD interaction and results in a deeper attractive potential (if chiral susceptibility is still 
negative) in chiral channel and a "larger" condensate! Such and such, the chiral condensate 
will become "larger" and "larger" as far as the chiral susceptibility become positive. In the 
limit m —>■ + , this process tends to be quasistatic, without change of free energy or entropy 
or energy, it means that a spontaneous quantum phase transition from Wigner-vacuum to 
Nambu-vacuum takes place! If the perturbative mass is negative, the situation is similar, 
but leads the Wigner-vacuum to another physically equivalent Nambu-vacuum which can be 
affirmed in the framework of D-S equation (in the chiral limit, we can always find two Nambu 
solutions for mass function with the same quantity and opposite sign, if the interaction in 
infrared region is strong enough [5], [11], [12] ). From Fig. 3, one can easily recognize that the 
susceptibility Fx is definitely negative if the chemical potential /i < 0.24 GeV, it becomes 
positive if the chemical potential fi > 0.24 GeV, and more significantly, the Fi is divergent 
and disconnected at fi = 0.24 GeV. Meanwhile Fig. 4 indicates that the susceptibility F2 
is positive if the chemical potential /i < 0.42 GeV and it changes to negative continuously 
as the chemical potential gets larger than 0.42 GeV. These behaviors manifest that the 
chiral symmetry breaking phase transits to the chiral symmetry preserving phase at chemical 
potential /1 = 0.24 GeV, and the phase transition is of first order, and that the diquark 
channel of Wigner-vacuum is stable while fi < 0.42 GeV, and the phase transition is of 
second order. Therefore we can reach a conclusion for the massless-QCD phase transition 
in the present truncated D-S equation approach as: there is a first order phase-transition 
of chiral restoration at a quite large chemical potential and color superconductivity phase 
emerges as a second order phase transition at a very large chemical potential where chiral 
symmetry has already been restored. 

To make our criterion more robust, we also investigate the coupling constant dependence 
of the chiral susceptibility F\ of Wigner solution at fi = and that of the Nambu dynamical 
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Figure 5: Coupling constant d dependence of the susceptibility F\ at zero chemical potential 
and zero momentum 



0.6 

> 
O 

— 0.5 



0.4- 



~ 0.3 



0.2- 



0.1 



0.0 



10 



20 



— I 

30 



parameter d [GeV" ] 



40 



I 

50 



Figure 6: Coupling constant d dependence of the dynamical mass function at zero chemical 
potential and zero momentum 



mass function M = Aim /An at n — 0. The obtained results are displayed in Fig. 5, 
Fig. 6, respectively. It is evident that if the interaction strength parameter d < 12.2 GeV -2 , 



10 



the chiral symmetry preserving phase (with dynamical mass M = 0) is stable. Only if 
the interaction in the infrared region is strong enough (in the present numerical case with 
uj = 0.4 GeV, the d should be larger than 12.2 GeV~ 2 ), can the chiral symmetry breaking 
take place. Meanwhile, it should be emphasized that while the critical coupling constant 
shown in Fig. 5 can be easily found by straightforward solving the D-S equation for dynamical 
mass function (see Fig. 6), the critical chemical potential shown in Fig.3 can not be obtained 
by solving the D-S equation for dynamical mass function (see Fig. 2). 

In summary, by solving the Dyson-Schwinger equations, we have studied the chemical 
potential dependence of the solutions of the D-S equation and of the Wigner-vacuum sus- 
ceptibility in the chiral and diquark channels as well as the effect of the interaction strength. 
It shows that if the chemical potential of the system is not very large and the interaction in 
the infrared region is strong enough, the QCD vacuum is in the chiral symmetry breaking 
phase. If the chemical potential gets larger and reaches a critical value, the chiral symmetry 
can be restored by means of a first order phase transition. If the chemical potential is much 
larger so as to arrive at another critical value, the color superconductivity phase emerges. 
Meanwhile the process of the chiral symmetry breaking is proposed to be a positive feedback 
with respect to the perturbation. Considering the criterions to characterize the chiral phase 
transition, we would like to mention that it is not advisable to introduce any approximate 
expression for free energy to judge which solution (Wigner solution or Nambu solution) is 
free energy favorable. First, no body knows how to construct an exact expression for free 
energy. It is more serious that an approximate expression can hardly satisfy the fact that 
Wigner-vacuum should always belong to the ground state Hillbert-subspace no matter the 
true vacuum is degenerate or not. Second, more importantly, it is not necessary to construct 
a free energy expression for understanding phase transition because the D-S equations (not 
solutions of D-S equations) contain naturally the full phase information of the dynamical 
system. For instance, as we have shown in this paper, the susceptibilities of the Wigner- 
vacuum can be widely used to understand dynamical symmetry breaking of the physical 
vacuum. 
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